Incoherent coincidence imaging and its applicability in X-ray diffraction 
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Entangled-photon coincidence imaging is a method to nonlocally image an object by transmitting 
a pair of entangled photons through the object and a reference optical system respectively. The 
image of the object can be extracted from the coincidence rate of these two photons. From a 
classical perspective, the image is proportional to the fourth-order correlation function of the wave 
field. Using classical statistical optics, we study a particular aspect of coincidence imaging with 
incoherent sources. As an application, we give a proposal to realize lensless Fourier-transform 
imaging, and discuss its applicability in X-ray diffraction. 
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Optical imaging techniques using classic light sources 
have been the primary tools for scientific research and 
industrial applications. In recent years, there has been 
an increasing interest in the field of quantum imaging, in 
which nonclassical states of light are used as light sources 
[ISSSHEHSSIIIlllJIIl Special attentions 
are focused on entangled-photon coincidence imaging 0, 

Humus. 

The role of entanglement in coincidence imaging leads 
to some debates now. The authors in Ref p stated that 
only quantum entangled sources can be used to realize 
coincidence imaging, and using classical light sources can 
not produce the image of the object. However, using 
classically correlated beams, the experiment performed 
in [l3j also produced a coincidence image. Moreover, in a 
recent preprint |l4|. it was shown using quantum theory 
that an object can be imaged via coincidence imaging 
with split incoherent thermal radiation. In this letter, 
we give a completely classical description of coincidence 
imaging and obtain a relationship between the intensity 
correlation in the detectors and the source. Especially, 
with a proper choice of the imaging geometry, we find it 
is possible to realize a kind of lensless Fourier-transform 
imaging by using an incoherent light source, which may 
be applicable for X-ray diffraction. 

An example of the setup of a coincidence imaging sys- 
tem is shown in Fig Q [y, ■ H the source S produces 
pairs of entangled photons, the produced photons are 
transmitted through a known (reference) optical system 
and an unknown optical system (test) which contains the 
object to be imaged. These two optical systems are char- 
acterized by their impulse response function h r (x,x r ) 
and h t (x,Xt) respectively. Two detectors D\ and D2 
record the intensity distribution of the test and refer- 
ence photons. The coincidence rate of photon pairs at 
these two detectors (G^ 2,2 -* (x r , Xt)) is proportional to the 
fourth-order correlation function of the optical fields, 

G {2 ' 2 \x r ,X t ) =<I{x r )I{Xt) >, (1) 

where < ... > means the ensemble average. In 




FIG. I: A setup of entangled-photon coincidence imaging. 
The source S emits pairs of entangled photons. One of the 
photons transmits through the test system ht(x,xt) which 
contains an unknown object, and the other photon transmits 
through a known reference system h r (x,x r ). Two detectors 
Dt and D r record the intensity distribution. The coincidence 
rate ,x t ) is measured to give a image of the object. 

entangled-photon imaging, the object can be ex- 
tracted from the marginal coincidence rate (I^ 2 \x r ) = 
Jdx t G^ 2,2 - 1 (x r , x t )) or the conditional coincidence rate 
(4 2) (x r ) = G< 2 ' 2 )(av,0)) @. Although the reference 
photons do not pass through the object, the object con- 
tained in the test system can be imaged at the reference 
detector. Such a nonlocal imaging technique may be use- 
ful for secure information transfer. 

Now, suppose the light source S is a classical light 
source, we use classically statistical optics to describe 
the coincidence imaging process. In the framework of 
fluctuating optical fields [15j, the forth-order correlation 
function G^ 2,2 ) (x r , x t ) relates to the optical fields in the 
reference and test detectors by 

G^' 2 \x r ,x t ) =< E* r (x r )E*{x t )E r (x r )E t (x t ) >, (2) 

where E r (x r ) is the optical field in the reference detector 
and Et(xt) is the optical field in the test detector. For 
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simplicity, we assume the source is quasi-monochromatic 
with a mean wavelength A. Only one transverse dimen- 
sion (x) is considered though the generalization to two 
transverse dimension is straightforward. 

If the optical field in the source is represented by E(x), 
the propagation of E{x) through two different optical 
systems leads to 



(3) 



Ek{x k ) = dxE(x)h k (x,x k ), 



where k = r,t. Note that h r , h t are deterministic func- 
tions, by substituting Eq. into Eq. (J2J, we have 

G { - 2 ' 2 \x r ,x t ) 

dx\dx' x dx 2 dx' 2 G^ 2 ' 2 '{x\, x'i, X2, x' 2 ) 
xh r (x 1 ,x r )h*(x' 1 ,x r )h t (x 2 ,x t )h*(x' 2 ,x t ). (4) 

where 

G ( - 2 > 2 \x 1 ,x' 1 ,x 2 ,x' 2 ) =< E*{x' l )E*{x' 2 )E{x 1 )E{x 2 ) > 

(5) 



is the fourth-order correlation function of the optical 
fields at the light source. 

Eq. |J2J| establishes the relation between the coincidence 
rate at the detectors and the correlation at the source. 
We need to know the properties of the light source to go 
further. In many cases, the field fluctuations of a clas- 
sical light source can be modeled by a complex circular 
Gaussian random process with zero mean [l5j, then 



G ( - 2 > 2 \x 1 ,x' 1 ,x 2 ,x' 2 ) 



G^ 1 \x 1 ,x' 1 )G^ 1 Hx 2 ,x' 2 ) + 
G^(x 1 ,x' 2 )G ( - 1 ' 1 Hx 2 ,x' 1 ),(6) 



where G^ 1 ' 1 )^, Xj) is the second-order correlation 
function of the fluctuating source field, represented 
by G^ 1 ' 1 ' (xi, Xj) =< E*(xi)E(xj) >, and satisfies 
G^ 1 )(x h x j )=[G^(x j ,x i )]*. 

Substituting Eq. @ into Eq. we get 



G^ 2 \x r ,x t ) 



dxidx' 1 G <,1 ' 1 \xi,x' 1 )h r (xi,x r )h*(x' 1 ,x r )) x ( dx2dx' 2 G < ' 1 ' 1 \x 2 ,x' 2 )ht(x 2 ,xt)hf(x 2 ,xt) 



+ (^j dx 1 dx' 2 G ( - 1 ' 1 \xi,x 2 )h r (xi,x r )hl(x 2 ,x t )^ X (^J dx 2 dx' 1 G ( - 1 ' 1 \x 2 ,x' 1 )h t (x 2 ,x t )h*(x / 1 ,x r 



< I r (x r ) >< It{x t ) > + 



dx\dx' 2 G^ 1 ' 1 ' {xi, x' 2 )h r (xi, x r )ht (x' 2l x t ) 



(7) 



The first term in the right side of Eq. J2J is the mul- 
tiplication of the intensity distribution at the reference 
and test detectors, and cannot be used to realize the co- 
incidence imaging [oj. |(j. However, if G^ 1 ' 1 ' (x%, x 2 ) is not 
factorable, the second term in the right side of Eq. has 
the similar form as in the entangled-photon coincidence 
imaging, apart from the presence of a phase conjugated 
/if instead of h t . Since a second-order correlation func- 
tion of a classical light source is factorable only when 
the source is fully coherent, we can perform the coinci- 
dence imaging using partially coherent or incoherent light 
sources. 

Let us introduce the intensity fluctuations in the two 
detectors 



AI k (xk) = h{xk)~ < Ik(xk) >, 



(8) 



in which k — r,t. The correlation between the intensity 
fluctuations at the reference and test detectors is 

< AI r (x r )AI t (x t ) > 



dx\dx' 2 G*- 1 ' 1 -* (xi, x' 2 )h r (xi 7 x r )h* t (x' 2l Xt) 



.(9) 



This correlation function is experimentalhy measurable. 
A similar result has been derived in Ref. [14|, but quan- 
tum theory is used in the derivation. 

Based on Eq. ©, we propose a scheme to realize lens- 
less Fourier-transform imaging by selecting proper h r and 
h t . 

Suppose the light source is fully spatially incoherent, 
then 



G^(x u x 2 ) =I(x 1 )S(x 1 -x 2 ), 



(10) 



where I(x) is the intensity distribution of the source and 
6(x) is the Dirac delta function. 

Further, the reference system contains nothing but 
free-space propagation from S to D r . Under the parax- 
ial approximation, the impulse response function of the 
reference system is 

h r (x,X r )= ... CXp <^ — -{X - X r ) 2 \ . (11) 



i\d r 
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where A is the source wavelength and k = 2tt/X is the 
wave number, d r is the distance between S and D r . 

The test system comprises an object at a distance di 
from S and a distance di from D t . The wave emitted 



from the source propagates freely to the object character- 
ized by the transmittance t(x'), then after transmission, 
it propagates freely to the test detector. The impulse 
response function of such a test system is 



, , s [,,e-' lkdl f-Mr. ,, 2 } . „e- lkd2 ( — vk , 

^ Xt) = J dx -Jxdv exp 1 m; {x ~ x ] ) t{x ) ^xdV exp \M^ {xt - x y 

Substituting Eqs. (|10I11I12|I into Eq. Q, after some calculations, we have 



< AI r (x r )AI t (xt) > 



dx'dx I(x)- 



z ikdr | in 2 | f 

-iXd r 1 Xd r 



lkdl f-*T, /xal , , N e- ifcd2 f-iTr. ,. 2 



(12) 



(13) 



If the source is large enough and the intensity distribution is uniform, we can regard I{x) — Iq, then Eq. (|13|l becomes 



< AI r (x r )AI t (x t ) >-- 



g—ik(di —d r ) 

dx' Jq • \ [ j — exp 



; '>< ,, 9 I , ,,e lkd2 { —in , 

" ' — ^''Ni» — i.n /i 2 
iACl2 I Act2 



■(a; r ~a;') 2 U(cc') 



iA(di — d r ) \ A(c?i — d r ) 

Selecting d±, c?2 and d r to satisfy d\ — d r = —d2, the quadratic terms of x' can be canceled, Eq. (|14(l has the form 



(14) 



< AI r (x r )AI t {x t ) > 



dx OT eXP lA^ ( *" r) J (x)eXP l A^2 



-'o 



T 



A^ 



r 



(15) 



where T(q) is the Fourier transformation of t(x'). So the 
correlation function between the intensity fluctuations at 
the reference and test detectors is the Fourier transfor- 
mation of the transmittance of the object. We note that 
the appearance of h\ rather than h t in Eq. (JJJ) allows this 
particular result to be obtained without the use of a lens 
anywhere in the system. 

If we measure the conditional correlation function of 
the intensity fluctuations by using a point-like test de- 
tector located at x t — 0, 



AI^ 2 \x r ) =< AI r {x r )AI t (0) >, 



(16) 



it will generate a image recorded in the reference detector 
but contains the information of the object. Equations 
(HSJ and JTBJ| then yield 



A^ 



T 



-2irx r 
Xd 2 



(17) 



So we found that, under the conditions of a large, uni- 
form, fully incoherent light source, without any optical 
instruments (such as lens) in the reference and test sys- 
tems, using a point-like detector D t and an array of pixel 
detectors D r , such a coincidence imaging system realizes 
the function of Fourier-transform imaging. 



Due to the success of the oversampling approach, co- 
herent X-ray diffraction imaging has attracted much at- 
tention recently [16-18]. However, several factors still 
limit the imaging quality. Because it is very difficult 
to fabricate optical components (such as lenses) that 
function in the X-ray regime, free-space propagation is 
used to obtain the diffraction pattern. Also, it is well 
known that currently used X-ray sources are generally 
incoherent. To achieve the spatial coherence needed to 
form high-quality diffraction patterns, such X-ray sources 
must be small and far from the object [19]. These require- 
ments decrease the illumination efficiency and necessitate 
the use of high brightness sources such as synchrotron 
sources. 

The lensless Fourier-transform imaging proposal given 
in this letter can overcome these difficulties. In fact, the 
image obtained in Eq. 117fl is exactly the diffraction in- 
tensity pattern of the object. Since there is no require- 
ment on the fully coherence, any kind of X-ray source 
can be used to realize X-ray diffraction imaging. As our 
method is insensitive to phase fluctuations of the source, 
the signal-to-noise ratio will be better than that achieved 
in direct diffraction imaging with an incoherent (or per- 
haps even partially coherent) source. So the incoher- 
ent coincidence imaging technique is applicable for X-ray 
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diffraction. 

Finally, we would like to discuss the effects of the time 
response of the detectors on our new imaging scheme. 
Generally, the intensity correlation < I r (x r )I t (x t ) > is 
not exactly measurable due to the finite time response of 
the detectors. Instead, we can only measure 

< l' r {x r , t)l' t {x t ,t + t) > 

rt+T R /2 i-t+r+T R /2 

= V / / < I r (x r ,t')I t {xt,t") > 

Jt-T R /2 Jt+T-T R /2 

xdt'dt". (18) 

where we write down the time dependence explicitly. In 
Eq. H18|l , r) is a coefficient and Tr is the average time re- 
sponse of the detectors. Since in our imaging scheme, the 
free space propagation distance in h r and ht are equal, 
the time delay r « 0. In X-ray range, Tr is much larger 
than the coherent time of the fluctuated fields, so the 
integration of Eq. (|18(l will be proportional to the equal- 
time intensity correlation < I r (x r )It(xt) > [13 ■ Actu- 
ally, in recent synchrotron radiation experiments, spatial 
intensity correlation have been measured by using slowly 
response detectors |2l|. The key point is that, since only 
spatial intensity correlation is concerned in our imaging 
scheme, using slowly response detectors will screen the 
temporal fluctuation and measure the spatial intensity 
correlation only. 

In conclusion, we have shown that a classically incoher- 
ent light source can be used to realize coincidence imag- 
ing based on the measurement of the correlation between 
the intensity fluctuations. Our treatments are fully clas- 
sical and do not use quantum theory. As an application, 
a scheme to realize lensless Fourier-transform imaging is 
described, which may be very useful in X-ray diffraction 
imaging. These results will be generalized to three di- 
mensional and the effects of source distribution or other 
imperfections will be considered in future works. 
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